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BAKER BOTTS, L.L.P. 
30 ROCKEFELLER PLAZA 
NEW YORK, NEW YORK 10112 



TO ALL WHOM IT MAY CONCERN: 

Be it known that I, Kenneth Shepard, a citizen of the United States of America, 
residing in City of Ossining, County of Weschester, whose post office address is 46 Snowden 
Avenue, Ossining, NY 10562, respectively, have invented 

METHODS FOR ESTIMATING THE BODY VOLTAGE OF DIGITAL 
PARTIALLY DEPLETED SILICON-ON-INSULATOR CIRCUITS 

of which the following is a 



The U.S. Government has certain rights in this invention pursuant to the terms of the 
National Science Foundation award number CCR-97-34216. 



The present invention relates to partially depleted silicon-on-insulator technology, 
and more particularly to techniques for estimating the body voltage of digital integrated circuits and 
partially depleted silicon-on-insulator ("PD-SOI") technology. 



SPECIFICATION 



NOTICE OF GOVERNMENT RIGHTS 



BACKGROUND OF THE INVENTION 



Field of the invention. 
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II. Description of the related art. 

Silicon-on-insulator ("SOI") technology has long been used for radiation-hardened 
or high- voltage integrated circuits. In more recent years, SOI has also been considered for us as a 
technology for high-performance, low-power deep-submicron digital integrated circuits. For digital 
5 applications, fully-depleted devices have been largely abandoned in favor of partially-depleted 
technology, because of the difficulty in controlling the threshold voltage of fully-depleted thin-film 
transistors. Partially-depleted SOI ("PD-SOI") technology has two main advantages for digital 
applications: the reduction of parasitic source-drain depletion capacitances and the reduction of the 
reverse-body effect in stack structures and pass-transistor logic. 

|© At the device and circuit level, however, the floating body effect in PD-SOI has posed 

and continues to pose major challenges in the successful use of this technology. As discussed in 

O Pong-Fei Lu et al., "Floating-body Effects in Partially Depleted SOI CMOS circuits, 32 IEEE J. 

fy Solid-State Circuits 1241 (1997), there is a parasitic bipolar effect which can result in noise failures 

if not correctly considered. In addition, there can be large "uncertainties" in the body potential, and 

15 consequently the threshold voltage, of devices due to unknown past switching activity. For many 

B ■- ■ 

yl circuits, the design margining required to protect against this uncertainty erodes all of the potential 

fj 

yjjj performance advantage under nominal operation. In addition, for many circuit styles in which noise 
margin is strongly determined by threshold voltage (e.g. dynamic circuits), considerable overdesign 
for noise can also result from conservative body-voltage margining. 

20 Several circuit design techniques, such as those described in D.H. Allen et al., "A 

0.20 ^m 1 .8 V SOI 550 MhZ 64b Power PC Microprocessor with Cu Interconnects," Digest Tech. 
Papers, ISSCC, pp. 438-39 (1999), and C.T. Chuang et al., "SOI Digital CMOS VLSI - A Design 
Perspective," 36th ACM/IEEE Design Automation Concerence, pp. 709-714 (1999), attempt to 
contain the noise impact of the parasitic bipolar current and the delay and noise-margin variation due 

25 to the floating body. Some of these design techniques, such as predischarging internal nodes of the 
nFET pulldown stack in domino logic to avoid parasitic bipolar currents, are quite counter to design 
practice in bulk. 
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In addition, previous circuit-level modeling work on PD-SOI has focused either on 
device issues, or on delay and noise effects due to the floating-body effect evident for particular 
circuits under periodic stimulus, such as pulse stretching, frequency-dependent delay time. See, for 
example, J. Gautier et al,. "On The Transient Operation of Partially Depleted SOI NMOSFET;s," 
16 IEE Electron Device Letters 498 (1995), and R. Puri et al, "Hysteresis Effect in Pass-Transistor- 
Based Partially-Depleted SOI CMOS Cirsuits," Proc. Int'l SOI Conf. 1 998. However, neither such 
modeling work nor the circuit design techniques discussed above are able to accurately estimate the 
upper and lower bounds of the body voltage of PD-SOI circuits in order to characterize the important 
electrical characteristics of the circuit. Accordingly, there exists a need in the field to accurately 
estimate the upper and lower bounds of the body voltage of PD-SOI circuits, taking into account the 
past history of the circuit, so that certain important electrical characteristics, such as delay and noise 
margins, of the circuit can be properly considered without requiring excessively long simulatrons. 

SUMMARY OF THE INVENTION 

An object of the present invention is to provide a technique for accurately accurately 
estimate the upper and lower bounds of the body voltage of PD-SOI circuits. 

A further object of the present invention is to provide a technique for determining the 
body voltage of PD-SOI circuits which accounts for the past history of the circuit. 

Yet a further object of the present invention is to provide techniques for accurately 
determining important electrical characteristics,such as delay and noise margins, of PD-SOI circuits. 

In order to achieve these objectives as well as others that will become apparent with 
reference to the following specification, the present invention provides a method for statically 
estimating a body voltage of one or more transistors which form digital partially depleted silicon-on- 
insulator circuit. In one technique, one or more device models are obtained and abstracted to 
generate one or more simplified electrical descriptions, each corresponding to one of the transistors. 
The device models may be n or p channel Field Effect Transistor models. As part of this abstraction, 
the circuit topology is then checked to generate one or more sets of accessible states, with each set 
corresponding to one of the transistors and being indicative of whether under any allowable 
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switching activity, the source, gate, or drain could be high or low. Next, sets of reference state body 
voltage minima and reference state body voltage maxima are determined for each of the transistors 
based on corresponding simplified electricaldescriptions and corresponding sets of accessible states. 
Finally, a target state body voltage minima and a target state body voltage maxima are ascertained, 
one for each transistor, based on the determined sets of reference state body voltage minima and 
reference state body voltage maxima. This target state reflects the appropriate initial condition for 
the electrical simulation in question. The electrical simulation in question also determines whether 
the minima or maxima is used for a given FET. 

Preferably, the abstracting step includes abstracting each device model to obtain 
displacement voltages d { , steady-state reference voltages Vj Zero , and forward biased reference 
voltages Vi forward , associated with each transistor state i for each such model. Using such abstracting, 
the determining step advantageously involves determining the minima and maxima reference state 
body voltage (V B ref ) min or (V B ref ) max for each transistor in the given circuit. With sucWa methodology, 
the present invention advantageously provides for ascertaining one or more target state body voltage 
minima, and one or more target state body voltage maxima, based on the determined reference state 
body voltage minima and maxima using either full uncertainty or accessibility analysis, and the 
displacement abstractions d r 

In an especially preferred technique, an additional step of performing active net 
tagging on each net in the circuit to determine whether any of nets will switch with regular frequency 
is performed prior to the determination of reference state body voltage minima and maxima. In this 

arrangement, the simplified electrical abstractions are thed/ V. zero and V forward combined with the 

results of active net tagging, and the determining step advantageously involves determining 
minimum and maximum reference state body voltages for each transistor in the circuit based on 
corresponding displacement voltages d i9 corresponding sets of accessible states. With this 
methodology, the present invention also calls for ascertaining one or more minimum target state 
body voltage minima and one or more target state body voltage maxima, based on the determined 
reference state body voltage minima and maxima using modified accessibility analysis. 
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In another preferred technique, the device models are also abstracted to obtain time 
constant characterizations, and an additional step of calculating signal probabilities and timing 
windows by Boolean and timing analysis is performed prior to the determination of reference state 
body voltage minima and maxima. In this arrangement, the determining step advantageously 
involves determining the sets of reference state body voltage minima and maxima based on 
corresponding displacement voltages d i5 the time constraints, corresponding sets of accessible states, 
and the results of the calculated signal probabilities and timing windows. With such a methodology, 
the present invention calls for ascertaining the target state body voltage minima and maxima, based 
on the determined reference state body voltage minima and maxima using probabilistic analysis. 

In order to take advantage of the ability to accurately determine target state body 
voltage minima and maxima, the present invention provides techniques for analyzing an electrical 
property of a digital partially depleted silicon-on-insulator circuit. In one preferred methodology, 
the technique involves ascertaining a target state body voltage minimum and a target state body 
voltage maximum for each transistor in the circuit, establishing an initial condition for the circuit 
by selecting either the target state body voltage minimum or the target state body voltage maximum 
for each transistor, applying a voltage to the circuit, and measuring the sought after electrical 
property of the circuit. 

In accordance with the present invention, the target state body voltage minima and 
maxima are ascertained alternatively from full uncertainty estimation, accessibility analysis, 
modified accessibility analysis, or probabilistic analysis. If the sought after electrical property is a 
switching delay, the measuring step involves measuring a delay between a switching input and a 
switching output. Such delay measurements can be used as constituent simulations in static timing 
analysis. If the sought after electrical property is noise, the measuring step involves measuring noise 
on one or more nets in the circuit such noise measurements can be used as constituent simulations 
in static noise analysis. In its ability to accurately determine target state body voltage minima and 
maxima, the present invention fulfills the need to accurately estimate the upper and lower bounds 
of the body voltage of PD-SOI circuits, taking into account the past history of the circuit, so that 
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certain important electrical characteristics, such as delay, induced noise, and noise margins, of the 
circuit can be properly considered. 

The accompanying drawings, which are incorporated and constitute part of this 
disclosure, illustrate a preferred embodiment of the invention and serve to explain the principles of 
the invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 is a state diagram for a partially depleted silicon-on-transistor pFET; 
Figs. 2a and 2b are graphs displacements for an example technology; 
Fig. 3 is an illustrative drawing showing the application of modified accessibility 
analysis to an exemplary circuit; 

Fig. 4 is an illustrative diagram showing an average cycle divided into a series of time 

slices; 

Figs. 5a-c are graphs plotting the relaxation times r t for an nFET with a 2.5-V supply 
as a function of local body voltage; 

Figs. 6a and 6b are graphs plotting reference body voltage against effective time; 

Fig. 7 is a flow diagram of the preferred process implemented by the present 
invention; . . . 

Fig. 8 is an illustrative diagram showing a chain of inverters; 

Figs. 9a and 9b are illustrative diagram showing an exemplary circuit; and 

Figs. 10a - 13c are graphs illustrating the results of a comparison between the 
techniques described herein and a simulation for the an exemplary circuit. 
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DESCRIPTION OF PREFERRED EMBODIMENTS 

In order to better understand the embodiments disclosed herein, the PD-SOI device 
physics which underlie the present invention will first be explained. It should be noted that the 
following description is by way of example and not by limitation, with certain embodiments being 
described in order to best explain the principles of the present invention. 

The body potential of a PD-SOI Field-Effect-Transistor ("FET") is determined by 
capacitive coupling of the body to the gate, source, and drain, by diode currents at the source-body 
and drain-body junctions and by impact ionization, currents produced by current flow through the 
device (referred to as the on-state impact ionization current). The impact ionization currents have 
a strong supply-voltage dependence, decreasing with decreasing supply voltages and going to zero 
below 1.1 -V supplies. 

It is convenient to distinguish "fast" and "slow" processes. Fast processes can 
change the body potential on time scales on the order of or less than the cycle time, while slow 
processes require time scales much longer than the cycle time (up to milliseconds) to affect the body 
voltage. There are two fast mechanisms at work: switching transitions on the gate, source, or drain 
which are capacitively coupled to the body (which will hereinafter be referred to as coupling 
displacements), and forward-bias diode currents across source-body and drain-body junctions with 
voltages exceeding the diode turn-on voltage (which will hereinafter be referred to as body 
discharge). The slow processes involve charging or discharging the body through reverse-biased 
or very weakly forward-biased diode junctions and through impact ionization. 

Figure 1 is a state diagram abstraction for an nFET which is useful to model the 
switching history determining the body voltage of a particular device. As a usually dynamic circuit 
node, the floating body has "memory." The states 1 10, 120, 130, 140, 151, 152 denoted with solid 
circles represent "static" states, states in which the FET can be stable. The states 161, 162 are 
"dynamic" states, which are only present transiently during switching events. For example, state 1 1 0 
corresponds to the case in which the gate is high and both the source and drain are low. Arrows 
indicate possible state transitions produced by switching events in the circuits containing these FETs. 
These switching events can represent transitions from the logic state at the end of the previous cycle 
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to the logic state at the end of the current cycle or can represent hazards that occur transiently within 
a cycle. States 151 and 152 can usually be treated equivalently as the same state 150; similarly 
states 161 and 162 can usually be treated equivalently as the same state 160. 

While Figure 1 applies to an nFET, is should be noted that the state diagram of a 
pFET is similar, in which the gate is high rather than low in states 130, 140, and 1 50; and low rather 
than high in states 110, 120, and 160. 

If the device is allowed to remain in one state for a very long time, the body voltage 
in each state will achieve a dc value, hereinafter referred to as si. The dc voltages in states 110 and 
130 (s, and s 3 ) are zero, while the dc voltages in stages 120 and 140 (s 2 and s 4 ) are given by the 
supply voltage. s 5 is determined by the steady-state balance between a weakly forward-biased 
junction drawing current from the body and a reverse-biased junction, leaking current to the body, 
enhanced by GIDL currents. Similarly, s 6 is determined by the steady-state balance between a 
weakly forward-biased junction drawing current from the body and charging curreht due to reverse, 
leakage of the other diode junction and on- state impact ionization. . i v : \ , 

As an example, commercially available BSIM3SOI models for an IBM partially- 
depleted SOI technology, having a 0.25 jum effective channel length, 5 - nm gate oxide, 350 - nm 
back oxide, and 140 - nm thin silicon film are considered herein. Those skilled in the art will 
understand that the present invention applies with equal force to other commercially available and 
custom fabricated partially-depleted SOI technology. 

Using exemplary supply voltages of - 1 .0V and 2.5 V, Table I lists the values of Sj 
for each of the six states. The remainder of the information listed in Table I will be explained below. 



NY02:266402.3 



Page -8- 



AP 32158-070050.1280 



TABLE I 
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1.77 


2.41 


0.0 


0.0 


-0.6 


0.0 


0.53 


1.13 


0 


0 


-0.6 


2 


2.5 


2.50 


3.1 


2.5 


0.73 


0.12 


1.0 


1.0 


1.6 


1.0 


0.469 


-0.127 


3 


0 


2.15 


2.85 


0.0 


-0.92 


-0.98 


0 


0.87 


1.5 


0 


-0.75 


-0.907 


4 


2.5 


3.43 


3.5 


2.5 


0.35 


-0.30 


1.0 


1.75 


1.906 


1.0 


0.14 


-0.5 


5 


0.437 


2.58 


2.79 


2.068 


-0.045 


-0.23 


0.278 


1.12 


1.47 


0.722 


-0.118 


-0.47 


6 


0.711 


1.60 


1.60 


1.976 


1.24 


1.24 


0.357 


0.61 


0.85 


0.722 


0.467 


0.17 



In the absence of body discharge, the coupling displacements that occur with each 
transition in Figure 1 are completely reversible on "fast" time scales; that is, if one begins in state 
110 and traverses the state diagram, returning to state 1 10 on a time scale much faster than any of 
the "slow" leakage mechanisms, the body voltage on return will be the same as the initial body 
voltage, a simple result of charge conservation. Because of this, the charge stored on the body can 
be represented as the value of the body voltage in one particular state of Figure 1 (hereinafter 
referred to as the "reference" state), which is chosen to be state 120 for the nFET and state 1 10 for 
the pFET. 

From this reference body voltage ), we can then determine the corresponding 
body voltage in each state i \y l B ) according to the relationship: 



V B= V B f+ V V B f ) (1) 

From ( 1 ), The displacements, 1 ^ B ' , are explicitly shown to be dependent on the reference body 
voltage because of the strong voltage-dependence of the capacitances of the source-body, drain- 
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body, and gate-body. d t is independent of device width W because the dominant components of 
capacitance scale proportionately with W. 

Referring next to Figures 2a and b, the displacements are shown as a function of 

yref 

b for the example technology at the 2.5-V supply, with Figure 2a illustrating an nFET and Figure 

5 2b illustrating a pFET. In many cases, these values are not easily determined with direct simulation 
of a transition from the reference state. For example, reference body voltages for the nFET greater 
than about 3.1 V would result in strongly forward-biased body- drain and body-source junctions 
making it difficult to distinguish the displacement due to the transition from state 1 20 from the body 
discharge. As a result, these curves are instead determined by isolated evaluation of the MOS C-V 
IQQ model. 

fi With the reference body voltage as a "state-independent" way of representing the 

% E charge trapped on the body, each state / in Figure 1 may be characterized by two values of this 

!>f reference voltage, V. zer and yf ormrd Referring again to Table I, y* ero represents the steady-state 

^ value of the reference body voltage, achieved by remaining in state i for a long time. This follows 

s 

1 S3 immediately from the s, in each state: 

m 

0 s^V^ + djiV*" 0 ) W 

O yf° rward represents the value of the reference body voltage for the nFET (pFET) to which the body 

would be very quickly pulled down (up) as a result of body discharge (charge), if state / were 

accessed with a higher (lower) reference body voltage than yf° nvard m The values shown in Table I 

20 presume that the fast body discharge will bring the forward-biased- junction bias down to a turn-on 
voltage of 0.6 V; it is also important to note that fast body discharge can trigger parasitic bipolar 
leakage between source and drain for FETs in state 150. This means, for example, that if a FET 

which reached a dc steady-state in state 140 (with a of 3.3 V) switches into state 120, the 
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reference body voltage will quickly discharge to yf° rward = 3.1 V. If the FET subsequently 
remains in state 120 for a long time, V$ ef will eventually decrease to V^™ = 2.5V. 

An important qualitative difference in the values between V. zero and v. fonmrd values 

the 2.5-V and 1.0-V cases is the 6 state. At a supply voltage of 2.5 V, v/ orward = V 6 zero (which 

is not the case at the 1 .0- V supply) because of the dominating effect of the on-state impact ionization 
current. This current is so large that a strongly forward-biased junction is required to balance it in 
steady-state. 

In accordance with the present invention, once the dj curves and the 
V. zer0 and v. fonvard values for each state are obtained, the target state body voltage minima and 

m 

10fl target state body voltage maxima are estimated. Four illustrative estimation techniques, full 

E s 

j^j uncertainty estimation, accessibility analysis, modified accessibility analysis, and probabilistic 
pjj analysis, may be used to ascertain approximate values of the target state body voltage minima and 

s target state body voltage maxima. Each is discussed in turn. 

Q 

m In the first technique, full-uncertainty estimation, no knowledge of the switching 

1 fy activity of the circuit is assumed. The maximum and minimum possible values of the body voltage 

□ that cover all possible stimulus and history are chosen. A state is accessible if the circuit topology 

D 

allows the state to be visited. As an example, for the nFET of an inverter, those states with the 
source high would not be accessible, because the source of the nFET is tied to ground. If A 
represents the set of all such accessible states, including possibly the dynamic state 160, the 
20 minimum and maximum body voltages are given by (3) and (4): 

r yzero (3) 

(V ref = max'; v } 

v B ' max jeA 



- . yzero 

(v re J \ = mm J (4) 

{y B *min jeA K J 
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For the second technique, accessibility analysis, it is assumed that the circuit is under 
steady switching activity such that every accessible state is visited with reasonable frequency, i.e., 
on a time scale that is faster than the "slow" body-voltage mechanism. Under such conditions, the 
y forward va i ues f or ^ e n FET (pFET) will cap the maximum (minimum) possible value of the body 

voltage, and values for the pFET will cap the minimum possible value of the body voltage. For the 
nFET, 

/yref^ = miny zero (5) 
v B ; min jeA j 

(V' ef ) «min(™x v?ero.min yforward (g) 

V ' max J eA J J eA static J V ' 



1(L 



while for the pFET, 

Ln 



nl 



(V r R ef ) . =max(^>f™;^ V{° rward ( 7 > 



{V™f) = max V™™ 



B 'max j (o\ 

^ j€A ' J \P) 

1 S The assumption here is that it does not matter how long an accessible state is visited; 



r 3 ™ 



jyj it will be long enough to discharge the body down to the forward-bias turn-on voltage of the source- 
y body or drain-body diodes. In general, state 1 60 is omitted from this "accessibility" analysis ( A static 

Q is used in the above equations) because state 160 is a switching state and, therefore, cannot be 

□ 

assured to meet this criterion. In those cases in which the relaxation in state 160 is extraordinarily 
20 "fast" because of a large on-state impact ionization current, state 160 can be safely included in the 
accessibility analysis. 

While accessibility analysis does not require detailed switching knowledge, it does 
require that there is enough switching activity that every accessible state is visited with a minimum 
frequency. At times, this, too, may be difficult to ensure. The third technique, which shall be 
25 referred to as modified accessibility analysis, permits the marking of signals in the design as "active" 
signals, i.e., signals that are assured to switch with regular frequency. As those skilled in the art will 
appreciate, clock nets can usually be tagged. This active net tagging establishes a set of constraints 
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that must be satisfied by a modified accessibility set of states, A ma . These are states that must be 
visited as a result of the "active" net switches. The body voltage is determined by equations 
identical to those used for ordinary accessibility analysis, except that A static is replaced by A ma . Thus, 
for the nFET, 



(V£pf) . = min V? er0 (9) 
j^A J 

(V ref ) = min (max V zer0 ) ; min V -forward (10) 
B max j j 

j e A j g A 

ma 

while for pFET. 

(V« f ) . = max(min V™ ro \ max v/ ormula ) nn 
B 'min v . j . j (11) 

/ e A / € A 

. . J ma 

(V^ ef ) = max V zero 

v B 'max / (lz) 

A ma will be the accessibility set for the nFET (pFET) which satisfies the constraints while ensuring 
the maximum (minimum) value of the body voltage. As in accessibility analysis, modified 
accessibility analysis explicitly excludes the dynamic state 1 60 from A ma . These constraints are 
derived from the following four rules: 

(1) If the gate of the FET is an active net, then there must be a state in A ma with 
the gate high and a state in A ma with the gate low. 

(2) If the source of the FET is connected to ground through a succession of 
transistors all of whose gates are active nets, then A ma must contain a state with the source low. If 
connected to supply, then A ma must contain a state with the source high. 

(3) If the drain of the FET is connected to ground through a succession of 
transistors all of whose gates are active nets, then A ma must contain a state with the drain low. If 
connected to supply, then A ma must contain a state with the source high. 
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(4) Same-signal correlations must be considered and can result in stronger 
constraints. For example, if gate of the current FET is marked as an active net, but must be high for 
a path to ground from the source to be present through a path of active-gate FETs, then A ma must 
contain a state in which both the gate is high and the source is low. 

Referring to Fig. 3, modified accessibility analysis is illustrated with reference to an 



exemplary circuit. Nets C and B (marked with the arrows) are active nets. For transistor Ml, 
therefore, there must be a state in A ma with the source and gate both high, with the gate low, and with 
the drain high. 



TABLE II 


nfet pfet 

T r zero jr forward . jrzero -rr forward 

i si v i si y. ys 
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1.78 


2.37 


0.0 


0.0 1 


-0.6 


2 


2.5 


2.50 


3.1 


2.5 


0.73.6 


0.12 


3 


0 


2.16 


2.79 


0 


-0.9315 


-0.9505 


4 


2.5 


3.43 


3.45 


2.5 


0.349 


-0.225 


5 


0.421 


2.54 


2.72 


2.079 


-0.026 


0.225 


6 


0.560 


1.3 


1.36 


2.078 


1.373 


1.15 



Table II lists values of V ( Zer0 and V- t forward for the nFET and pFET in the circuit of Fig. 
3. To find A ma that satisfies these constraints and achieves the maximum value of , refer to 

Table II and begin with the state of largest v/ orward . Adding state 1 40 to A ma satisfies the drain-high 

constraint and the gate-low constraint. State 120 satisfies the constraint that there must be a state 

with the source and drain both high. Therefore, A ma = {140, 120} and ( V B f ) m2LX = 3.1 V. 

The fourth and most complex technique presented herein for estimating the target 
state body voltage minima and target state body voltage maxima shall be referred to as probabilistic 



NY02:266402.3 



Page -14- 



AP 32158 - 070050.1280 

analysis. In this technique, the estimates provided by the simple uncertainty modes are tightened 
with stochastic techniques and timing information from static timing analysis. Thus, the behavior 
of the body over a long period of time is determined by reference to an average cycle repeated over 
and over. 

5 In probabilistic analysis, two average cycles are characterized, one to minimize and 

one to maximize the body voltage. Referring to Figure 4, such an average cycle is shown, divided 

into a series of time slices t*f , which characterize the amount of time per cycle on the average that 

the FET spends in state /. The time slices sum to the cycle time, in accordance with equation 13: 

6 e ff (13) 
10 i cycle 

O 
Hp 

m 

if1 1 ; In each time slice /, the reference body voltage is relaxing to the target value 

j_ y zero w i t k ^ characteristic relaxation time denoted by v r The body-voltage-dependent time 

'-4 ■ 

q constants for this relaxation are denoted by the ^ and can be captured from the device models as part 
1^ of a technology precharacterizationm, which may be a piecewise-linear representation of the 



!i' logarithm of r- t as a function of V™f , since r } varies over several orders of magnitude. Moreover, 



y 
u 

J3 

R Tj is independent of device width W because both the body currents and body capacitance scale 



w proportionately with W. 

Except for the fast discharge associated with a source-body or drain-body junction 
20 that becomes strongly forward biased, these time constants are much larger than t cycle , and any 
voltage change during a single time slice would be imperceptible in Figure 4. From the foregoing, 
the reference body voltage at the end of the cycle F„ +/ can be related to the reference body voltage 
at the beginning of the cycle V n by equation 14: 

eff/t -- 6 ~m -t.eff/t- , eff/t. (14) 
i=l 



25 V n + r V <4r'i VI, V r°<'-e 1 , ) e I^ i+| t j 
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The steady-state solution of this difference equation (in the approximation that the z; are much 
greater than t cyde ) is given by equation 15: 

If. i (t e ff ITi)V zero 
yref l ~ ] K * L? (15) 

B Zf.i tfff/Ti 



Because of the y™f -dependence of the x i? these equations must be solved self-consistently. 

Figures 5a-c plot the relaxation times r t for the nFET at the 2.5-V supply as a 
function of the "local" body voltage (that is, the body voltage in state / rather than the reference body 

lQL voltage). The relaxation times x, and x 3 , for states 110 and 130, peak around a zero body voltage 

B 

<U and then decrease for negative body voltage as the source-body and drain-body leakage currents 

fff 

J3 increase. The decrease for positive body voltage is even more substantial as the junctions become 
forward biased and draw significant current. The same trends are observed for the relaxation times 

S! x 2 and x 4 , for states 120 and 140, where the relaxation times peak around the supply voltage, 

\% indicating that the source and drain are both at the supply voltage for states 120 and 140. 

O 

^ The relaxation time x 5 for state 1 50 peaks around zero. For voltages above the peak, 

y » 

D x 5 decreases strongly as the body-source junction becomes forward biased. For voltages below the 

. E~E 

p peak, x 5 also decreases as both the body-source and body-drain junctions become more reverse- 
™ biased. x 5 peaks at a much smaller time-constant value than Xj through x 4 because even at a body 
20 voltage near the peak, one of the junctions (the drain-body junction) is strongly reversed biased and, 
therefore, leaking considerable current. 

For state 160, the time constants are considerably smaller (factor of 1000) than they 
are in state 150. This is due to the contribution of the on-state impact ionization current to the 
leakage currents from the floating body; a very low relaxation time will mean that devices that have 
25 any time in state 1 60 in their switching transients will have their body voltages effectively "pinned" 
at s 6 . At 1 .0 V, by contrast, x 5 and x 6 are comparable in magnitude. 



NY02:266402.3 



Page -16- 



AP 32158 -070050.1280 

In probabilistic analysis, the maximum and minimum possible values from 
Equation 1 5 are ascertained upon consideration of all of the allowable t e f . If the Tj in Equation 
15 were constant and equal, then to maximize one would simply maximize the time in those 
states with higher y. zero (subject to constraints on the maximum and minimum possible values of 
the ff imposed by the stochastic analysis); similarly, to minimize one must maximize the time 
in those states with lower y. zero . The x { however, are not constant, but because they have the 
property that they decrease rapidly for V™ f different from y zero states with the largest 

continue to have the biggest "pull" and the y zero "ordering rule" continues to hold. 

To formalize this, the rank of the state / is defined as an integer ^indicating the priority of states to 
maximize the body voltage. Higher ri means higher priority. For example, from Table I, for the 
nFET at 2.5 V, r 4 = 6, r 5 = 5, r 2 = 4, r 3 = 3, r } = 2, and r ( = 1. A complementary rank, r, , can also be 

defined for minimizing ; these priorities are simply the reverse of the maximum case. Such 

ordering, which favors states with high source or drain voltages to maximize the body voltages while 
favoring states with low source or drain voltages to minimize the body voltages, is technology- 
independent except for the relative positions of state 150 and 1 60, which depend on both technology 
and supply voltage. 

Referring next to Figures 6a-b, cycle times for the nFET at the 1 .0-V supply is shown 
as an example of how the steady-state body voltage depends on the tf f . The cycle time t cyde is 10 
nsec. Figure 6(a), assumes a hypothetical case in which only states 1 10 and 140 are accessible, and 
plots the reference body voltage as a function of t*f given that t e f = t eyC k - t e f . 



V ref _ V; 
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The values of V zero ; V/ or w ; V 4 zero are also noted. Given the accessibility of states 

110 and 140, the range from V zer0 to V x fonmrd (denoted with the arrow), would be given by 
accessibility body-voltage estimation. The probabilistic body voltage estimate becomes 
V zero as t*f -> 0 and monotonically increases with increasing t e f (r 4 > n). As -> tc y de(md 

5 t e f -» 0). rapidly increases to V™ ro . can only increase above v/ orward , with ff =0 

on the scale of the "fast" forward-biased discharge of the body. 

In Figure 6(b), the case in which only states 1 1 0 and 120 are accessible, is illustrated. 

Vg f increases monotonically as increases {r 2 > rj. In this case, F 2 2ero < yf° rward so that 

3 £ 

~ V^as t!f -> 0 does not actually reach the upper bound of accessibility analysis. 
10^ To determine allowable values of tf f for probabilistic analysis, information about 

fy 

S! the circuit environment of each transistor in the circuit, both logical and temporal, is required. The 

O - 

B logical environment of each FET is therefore characterized by a set of signal probabilities which 

o 

rf determine the possible states of the source, gate, and drain of the transistor at the end of a cycle. For 

O the nFET, it is preferred that seven signal probabilities are used, as follows: 

1|3 (1) P(G) = probability that at the end of a cycle the gate is driven high; 

0 _ 

(2) P(Z)|G) = probability that at the end of a cycle the drain is driven high 

given that the gate is low; 

(3) P( D^G) = probability that at the end of a cycle the drain is driven low given 
that the gate is low; 

20 (4) P(S G) = probability that at the end of a cycle the source is driven high 

given that the gate is low; 
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(5) P(S G) - probability that at the end of a cycle the source is driven low 

given that the gate is low; 

(6) P(S|G) = P(D|G) = probability that at the end of a cycle the drain (and source) 
are driven high given that the gate is high; and 

5 (7) P(S\G) = P(D\G) = probability that at the end of a cycle that the drain 

(and source) are driven low given that the gate is high. 

The corresponding probabilities for the pFET are: P(G), P(S|G), P(S\G), P(D|G), 

P(D\G) 9 P(S|G) = P(D\G), P(s\G)= P(d\G). The term "driven low" means that there 

: s 

S is a path to ground; "driven high" means that there is a path to supply. The following conditions on 

lQO the conditional signal probabilities must hold: . _ 

m P(D|G) + P(D|G) <1 

tj . . P(S|G) + P(SJG) < 1 

Q P(D|G) + P(D|G) <1 (16) 

P(S|G) + P(S|G) < 1 

IS 

p None of these have to sum to precisely 1 because of the possibility of floating nodes. 

~ With these signal probabilities, one can view the state diagram in Figure 1 as representing a Markov 
B process with six-by-six transition matrix A. Assigning different transition matrices to the minimum 
and maximum cases yields difference equations 1 7 and 1 8, where the six-by-six matrices 
20 A max and A min are given 1 9 and 20: 



/pmax ) _ a /pmax \ 



(18) 

K+\ 



(P mm ) K ^ = A • ( pmm ) 
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(19) 



, /rt max „ max „ max „ max n max n max > ~ 

' = ( p 1 P 2 P 3 P 4 P 5a P 5b > T 
max xfllllllj 



. /n min n min n min n min n min n min , T 

* min = ( p 1 P 2 P 3 P 4 P 5a P 5 ) T 

min xriiini; (20) 



10 In equations 17 : 20, (P,)£ is the probability of being in state i at the end of cycle £, 

~j and /?, is the probability of making a transition to state i and follows directly from the source, gate, 

Ul max _, - 

y3 and drain signal probabilities. For example, for the nFET, p 3 =(1 - P(G))P(S\G)P(D\G) 

M* 

n s 

j? while p = (1 - P(G))(l - P(S G))(l - P(D G)). The maximum case assumes that the 

« 

P floating node condition on the source or drain takes a high voltage value; the minimum case assumes 

y 5 

lf3 that the floating node condition on the source or drain takes a low voltage value. 

■ y3 

n Diagonalizing the matrix A and finding the eigenvector associated with eigenvalue 

W 1 (normalized so that the sum of the elements of the vector is 1) gives the steady-state values of the 
P h From these probabilities, one can calculate a set of thirty-six transition probabilities (that is, the 
probability that at given FET in a given cycle is transitioning from state i to state j) for both the 
20 minimum or maximum cases, in accordance with equation 2 1 : 

max max max min min min 

P = P P and P = P P . . (21) 
' -> j i J * -> J ' J v 7 

Next, the fraction of the cycle time ( that can be spent in each state / is 
determined as part of the transition j-k to maximize (or minimize) the body voltage among the set 
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of possible waveforms. To do this, early and late arrival times (rising and falling) for the source, 
gate, and drain of the FET under consideration (the target FET) are required. The seven early times 
are and are denoted as follows: 

(1) Srbe* • Earliest time the source of the FET can be driven high; 

5 (2) Sffi . Earliest time the source of the FET can be driven low; 

(3) Drise* * Earliest time the drain of the FET can be driven high; 

(4) D e j Q f . Earliest time the drain of the FET can be driven low; 

(5) G e f °il y . Earliest time the gate of the FET can be driven low; 

jjj (6) G™ r J y . Earliest time the gate of the FET can be driven high; and 

(7) G™ e rly . Earliest time the gate of the FET can be driven high. 

SJ- 

£3 There are comparable arrival times associated with the late state: 

b S%' e , Sff, , D'Z , D% , G% , and G% e e . For the source and drain arrival time, the target FET 

5 

P is assumed to be off. 

O 



Each transition has associated with it a set of arrival times necessary to make that 
1§^ transition, a transition set, denoted as T^ y For example, for the 1 10-* 120 transition for the nFET, 
the associated arrival time set is T^ 2 = {S ris& D ris J. For 1 10- 140 for the nFET, the associated arrival 
time set is T^ 4 ={G fa{b S ris& D ris J. The max and min operators which act on the transition sets are 
defined such that max(T i-j ) returns the largest of the early arrival times in the transition set, while 
min(Ti_j) returns the smallest of the late arrival times in the transition set. 
20 To indicate the states of Figure 1 involved in a cycle and to handle the possibility of 

hazards, the waveform in a cycle (in this case involving a transition from / to j) is denoted using a 
transition notation as follows: 

max T min 
/ > k > J 
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In this cycle, a hazard to state k occurs as part of the transition. The transition notation must involve 
only static states and indicates the amount of time spent in each of these static states as part of the 
transition. Specifically for this example, 

i -» j i -> j i -> j 

t . = max(T ); t , = min(T .)- max(T t ) ; and t = t , - min(T .). 
i i^k k k^J i^k J cycle k -» j 

Hazards are introduced when they act to increase (in the case that we are seeking the maximum body 
voltage) or decrease (in the case that we are seeking the minimum body voltage) the steady-state 
body voltage that would result from the particular waveform being repeated indefinitely. 

This transition analysis determines the t. J ~* k for the static states /. t^ k depends, 

in principle, on the number of switching events in a cycle which occur as a result of current flow 
through the target device. Each such switching event contributes an amount t switch to t . 7 ~* * where 

t switch is the approximate switching time of the FET. For the purposes of probabilistic body- voltage 
estimation, it is assumed that if state 1 60 is accessible, exactly one switching event occurs per cycle. 
In practice, the detailed results are not very sensitive to the number of assumed switching events or 
the exact value of t switch .. For cases in which the on- state impact ionization current is large, any 

• k produces pinning at V$ ef = y* ero yf° rward m p or cases i n which the on-state impact 
' 6 >0 

ionization does not dominate the current to the body in state 1 60, t*f is usually much smaller than 
the associated with the static states. As a result, probabilistic body voltage has very little 
sensitivity to the exact value of t*f . 
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From the tf~* k determined above, an effective amount of time (f 5^ ) on the average 

per cycle that the FET is in the state / is determined by two sets of t e f values, one to maximize and 
one to minimize the body voltage in accordance with equations 23 and 24: 

neff \max _ y pmaxft J"* ^Smax 
5 U i } j,k j^k^i > (23) 

/ t eff\min_y pmin / t ^\min 

l i ' j,k j->k^i > (24) 

Referring next to Fig. 7, the steps executed by a computer (not shown) executing 

1 (£3 software to implement the techniques described above will be described. In step 7 1 0, one or more 

yf| device models (n or p channel Field Effect Transistor models) are obtained. In step 720, the device 

models are abstracted to generate one or more simplified electrical descriptions, each corresponding 

to one of the transistors. As discussed above, the abstracting step may include abstracting each 
■ Nl ■ 
□ device model to obtain displacement voltages d { 721 , steady-state reference voltages Vi Zero 722, and 

1 forward reference voltages Vi forward 723, for each such model. Time constant characterizations 724 

JfJ are conditionally abstracted in the event that probabilistic determination is employed. 

a In step 730, the circuit topology is checked to generate one or more sets of accessible 

□ 

p states, with each set corresponding to one of the transistors and being indicative of whether a 
connection between a source or a drain of a transistor and either a power supply or ground exists. 

20 The conditional step 735 of performing active net tagging on each of net in the circuit to determine 
whether any of nets will switch with regular frequency is performed where modified accessibility 
analysis is employed. An alternative conditional step 736 of calculating signal probabilities and 
timing windows from the time constant characterizations by boolean and timing analysis is 
performed where probabilistic determination is employed. 

25 In step 740, sets of reference state body voltage minima and reference state body 

voltage maxima are determined for each of the transistors based on corresponding simplified 
electrical descriptions and corresponding sets of accessible states. As discussed above, this step may 
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include determining the sets of reference state body voltage minima and maxima based on 
corresponding displacement voltages d; and on corresponding set of accessible states 741, on 
corresponding displacement voltages d i3 corresponding sets of accessible states, and from the results 
of the active net tagging 742, or on corresponding displacement voltages d i? corresponding sets of 
5 accessible states, and from the results of the calculated signal probabilities and timing windows 743 , 
depending on the type of analysis being performed. 

Finally, in step 750, the target state body voltage minima and a target state body 
voltage maxima are ascertained, one for each transistor, based on the determined sets of reference 
state body voltage minima and reference state body voltage maxima. Full uncertainty estimation 
10 751, accessibility analysis 752, modified accessibility analysis 753 or probabilistic analysis 754 
O may be employed. 

IH In order to take advantage of the ability to accurately determine target state body 

jr voltage minima and maxima, techniques for analyzing an electrical property of a digital partially 
Ty depleted silicon-on-insulator circuit are now described. In one illustrative embodiment^ the 

SA 

1 53 technique involves ascertaining a target state body voltage minimum and a target state body voltage 
™ maximum for each transistor in the circuit, establishing an initial condition for the circuit by 

«| selecting either the target state body voltage minimum or the target state body voltage minimum for 

u 

y|3 each transistor, applying a voltage to the circuit, and measuring the sought after electrical property 

O 

p of the circuit. Techniques for analyzing two exemplary electrical properties, switching delay and 

20 noise, are now discussed. 

First, the body-voltage estimates described above are applied to a static transistor- 
level timing analysis engine, or a SOI Static Timing Analyzer ("SOI-STA"), in order to analyze the 
electrical property of switching delay. The results for full-uncertainty , accessibility and probabilistic 
analysis are considered, although modified accessibility could also be easily applied in the context 

25 of static timing analysis. The design is partitioned into channel-connected components (CCCs) for 
analysis, as is traditionally done in static transistor-level tools. SOI-STA utilizes a breadth-first 
search (BFS) of the resulting timing graph, which ensures that the arrival times on all of the inputs 
are known when the delay through the CCC must be calculated. This enables detailed body-voltage 
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estimates on all of the FETs of the CCC to be used in delay simulation. Delay propagation through 
a CCC occurs as a result of a single switching event on an input (i. e., simultaneously switching 
inputs are not considered). Depending on the degree of knowledge of past switching activity (no 
knowledge, quiescent, steady switching, steady switching with known signal probabilitiesand arrival 
times), one of the techniques discussed above (full-uncertainty, initial-condition, accessibility, 
detailed) is employed to determine the minimum and maximum possible value the reference body 
voltage can have for each FET of the CCC under analysis. For the given sensitization of the CCC 
for delay calculation, each FET is in a known state, the minimum or maximum body voltage of 
which can be determined by a displacement from the reference voltage. The body voltage values 
are then used as the "initial conditions" for the required delay simulation. 

Early-mode calculation for rising transitions sensitizes the nFETs of the pull-up path 
to be maximum and (to reduce the "fight" during switching) the nFETs of the pull-down path to be 
minimum; similar calculation for the pFETs of the pull-up path to be minimum and the pFETs of 
the pull-down path to be maximum. f ; , 

This same sensitization applies to late-mode fall transitions. Early-mode calculation 
for falling transitions sensitizes the nFETs of the pull-down path to be maximum and the nFETs of 
the pull-dup path to be minimum, and falling transitions sensitizes the pFETs of the pull-down path 
to be minimum and the pFETs of the pull-dup path to be maximum. 

The same sensitization applies to late-mode rising transitions. The fanout CCCs are 
included in the delay calculation to improve the delay accuracy. To maximize device capacitance 
for late-mode, the maximum nFET and minimum pFET body voltages are used. To minimize device 
capacitance for early-mode, the minimum nFET and maximum pFET body voltages are used. 

SOI-STA propagates full piecewise-linear ("PWL") waveforms. One of the 
complexities of BFS timing analysis is determining the late- and early-mode waveforms at the output 
of each CCC. It can sometimes be the case that the waveform with the maximum (minimum) delay 
(as measured by the 50-percent point) is not the waveform with the slowest (fastest) slew (as 
measured by the 20-to-80-percent rise-fall time). In these cases, a "hybrid" waveform is propagated, 
chosen with the largest slew as the late-mode waveform and translate it in time so that is has the 
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largest delay, or with the smallest slew as the early-mode waveform and translate it in time so that 
is has the smallest delay. SOI-STA also propagates signal probabilities using assumptions of spatial 
and temporal independence, borrowing from similar techniques in static power analysis. If 
probabilistic body-voltage estimation is used, once the signal probabilities and arrival times are 
known at the inputs of a CCC, these probabilities are translated into FET signal probabilities and 
arrival time values. 

Referring next to Figure 8, a 4-2 compressor circuit from a tree multiplier design is 
examined as an example of the advantageous use of the static transistor-level timing analysis 
presented herein, with the long delay path, I4-Ap-C-D-E-F-S, and the short delay path, I2-C-E-S, 
being considered. In the dynamic analysis, input waveforms of "50% duty cycle" are sued to 
sensitize these paths, corresponding to signal probabilities of 0.5 on all of the inputs in SOI-STA 
detailed steady-state analysis. In Table III, the initial-condition delays determined by circuit 
simulation are compared with the SOI-STA-determined full-uncertainty delays. 



TABLE III 




SOI-STA full-uncertainty 




SPICE initial 








Short Path 


Long Path 


Short Path 


Long Path 




rise 


fall 


rise 


fall 


rise 


fall 


rise 


fall 


12 


100 


100 


100 


100 


100 


100 


100 


100 


14 


100 


100 


100 


100 


100 


100 


100 


100 


Ap 






302 


491 






235 


315 


C 


301 


317 


866 


712 


326 


379 


476 


432 


D 






1057 


1079 






723 


553 


E 


526 


533 


1421 


1359 


1390 


762 


785 


943 


F 






1620 


1697 






1112 


936 


S 


672 


685 


1926 


1829 


1790 


837 


1104 


1174 
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Table IV does a similar comparison between the steady-state delays determined by 
circuit simulation (after more than 50,000 cycles of simulation) and the SOI-STA-determined 
probabilistic ("detailed")and accessibility steady-state delays. In all cases, the SOI-STA-determined 
delays bound the SPICE delays. One should also notice the considerable reduction in uncertainty 
between the full-uncertainty and probabilistic steady-state delays, as the component of this 
uncertainty due to body voltage variation is noticeably reduced. 



TABLE IV 

SOI-STA probabilistic SOI-STA accessibility SPICE steady-state 



Short Path Long Path Short Path Long Path Short Path Long Path 



12 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 , 


100 


100 


14 


100 


100 


100 


100 


100 


100 


100 


100 


100 i 


100 [ 


MOO ; 

i 


100 | 


Ap 






267 


409 






296 


449 






266 i 


362 i 


C 


464 


485 


706 


605 


389 


365 


804 


673 


488- ! 


:503 j 


594 , 


576 ; 


D 






897 


877 






1007 


999 






784 


794 


E 


970 


985 


1168 


1146 


733 


746 


1332 


1287 


1040 


1017 


1059 


1099 


F 






1363 


1409 






1540 


1590 






1246 


1352 


S 


1372 


1347 


1610 


1534 


1026 


986 


1812 


1737 


1384 


1365 


1508 


1483 



The body-voltage estimates described above are next applied to transistor-level static 
noise analysis, in order to analyze the noise immunity of the circuit. In this case modified 
accessibility analysis is considered for estimatics, although full uncertainty, accessibility, or 
probabilistic estimation could be applied to the circuit. As in the static time case, the approach 
involves decomposing the design into a collection of channel-connected components ("CCCs"), 
transistors that are connected together through their sources and drains. The maximum noise that 
is possible on each net is calculated as a time-domain waveshape. This worst-case noise considers 
all possible noise sources: leakage, charge-sharing noise, coupling though the interconnect, and 
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power-supply noise. This is done with a careful choice of vectors on the driving CCCs, referred to 
as the sensitization, which produces this worst case noise. 

Noise can also propagate from CCC-input to CCC-output (propagated noise). Noise 
failures are determined by the noise stability, a type of AC noise margin analysis, of each CCC given 
the worst case noise appearing at its inputs. This involves calculating the transient sensitivity of the 
output noise with respect to the dc-level of the input noise. While the special considerations 
associated with applying static noise analysis to PD-SOI circuits is discussed herein, the article by 
K.L. Shepard et al., "Harmony: Static Noise Analysis for Deep-Submicron Digital Integrated 
Circuit," IEEE Trans. CAD 1 1 32-1 1 50 (1 999), the disclosure of which is incorporated by reference 
herein, provides more details on static noise analysis generally. 

There are two important considerations for static noise analysis for PD-SOI circuits. 
The first is that all of the body voltages must be initialized as part of the CCC analysis used to 
. calculate/the. noise on each node as well as the :CCC analysis, used to determine the -noise stability 
of the gate. These body voltages are calculated' using the .modified accessibility analysis 1 described 
above The approach here is straightforward and depends on the noise type being calculated. Two 
types of noise on each CCC output, V H noise, which is noise that pulls the output down from the 
supply level, and V L noise, which is noise that pulls the output up from ground, are calculated. 

When calculating V H noise on a CCC output (or when verifying noise stability in the 
case that V H noise is introduced at the output), all of the devices in the pull-down stack are initialized 
to maximize the device strength by minimizing threshold voltages (maximum body voltages for 
nFETs and minimum body voltages for pFETs). This increases the strength of these devices in 
introducing noise. Similarly, all the devices in the pull-up paths are initialized to minimize device 
strengths by maximizing threshold voltages (maximum body voltages for pFETs and minimum body 
voltages for pFETs). This reduces the strength of these devices in maintaining the output at the logic 
high level. For V L noise, the situation is just the opposite. Devices in the pull-down path are 
weakened (minimum body voltages for nFETs and maximum body voltages for pFETs), while 
devices in the pull-up path are strengthened (maximum body voltages for nFETs and minimum body 
voltages for pFETs). 
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The second important consideration in the static noise analysis of PD-SOI circuits 
is that in calculating both propagated noise and in verifying noise stability, a parasitic bipolar 
leakage current may also be activated. Because of this, to achieve the worst-case, the arrival time 
of the input noise must be varied to maximize the output noise; that is, make sure that the noise 
5 introduced by parasitic bipolar leakage is superimposing maximally with the propagated noise 

Referring next to Figs. 9a and 9b, a simple circuit is examined as an example of the 
advantageous use of the noise analysis presented herein. Figure 9a illustrates a dynamic simulation 
of a circuit. Figure 9b shows the CCC decomposition used internally by the static noise analysis 
engine on the same circuit. In this partitioned circuit, the dynamic node is referred to as dyns and 
1 0 the output of the inverter as outs. 
O The static noise analysis results and dynamic simulation of the same circuit are 

ifp compared by vectors which correspond to the * " worst-case" behavior detected in the static analysis: • m - 
^ In the circuit of Figure 4(a), the coupling noise introduced by the switching of an aggressor %&^M^ * 
PJ or I) capacitiyely coupled to inputs A, C, and E, respectively, after a long period of switching histpty ! 
1 5^ is considered, called the prelude. Three different preludes are considered as shown in Table ^ beldwru* ' : 





TABLE V 


LB s 

o 


Prelude 






Node disposition 






P 




A 


B 


C 


D 


E 


F 


o 


1 


0 


0 


0 


0 


active 


0 




2 


0 


0 


active 


0 


0 


0 




3 


0 


0 


0 


0 


active 


active 



25 The clock period is lOnsec and the prelude is simulated for lOOmusec (so that the 

body has reached steady state) before the noise is applied. For the purposes of the dynamic 
simulation, a net that is noted as active in Table V is given a 0 - 1 transition at the middle of the 
evaluate cycle, every other evaluate cycle during the prelude. After the prelude, which is different 
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depending on the case considered, the "noise condition" shown in Table VI is applied, introducing 
V L coupling noise on node A for case 1, node B for case 2, and node C for cases 3 and 4: 



TABLE VI 


Case 


Prelude 


A 


B 


C 


D 


E 


F 


G 




H 


I 


1 


1 


0 


1 


0 


0 


0 


0 


fall 


1 


1 


2 


2 


0 


0 


0 


1 


0 


0 


1 


fall 


1 


3 


1 


0 


0 


0 


0 


0 


1 


1 


1 


fall 


4 


3 


0 


0 


0 


0 


0 


1 


1 


1 


fall 



This "noise condition" is identical to the one established by static noise analysis to produce the 
worst-case V H noise on node dynd in each case. 

The results of the dynamic circuit simulations are shown in Figures 1 0a through ;1 3a. 
In each of these graphs, the dynamic node dynd, the dynamic gate output node outs, the noise 
being injected on the gate, and the body voltage of the FET receiving the ; Y% noise on its gate is 
shown. The dynamic gate is given a " ' full-latch" keeper so that the failure of the case is immediately 
evident with the switching of this latch. Cases 1 and 3 show noise failures (as evidenced by the false 
switching of the latch), while cases 2 and 4 do not. In this dynamic case, the full prelude is 
simulated before the noise condition is applied. 

Figure 10b through 13b show the noise on node dyns computed by static noise 
analysis. It should be noted that in the static analysis, the feedback path is broken and replaced by 
a tie. In this case, these is no "prelude simulation." Instead, the body voltages of each FET are 
initialized by the tool with the maximum or minimum value as determined by modified accessibility 
analysis. The inputs noted as active in Table V (in addition to the clock) are considered as active 
nets in this analysis. In Case 1 , a significant parasitic bipolar current is also present; as a result, the 
tool adjusts the arrival time of the noise on A to produce the maximum peak noise at dyns. 

Figures 1 0c through 13c show the sensitivity of the domino-gate inverter output with 
respect to the noise calculated on dyns. Failures (as determined by sensitivity exceeding one in 
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magnitude) are correctly predicted for Cases 1 and 3. The noise is so large in Case 1 that the 
sensitivity has two " peaks" as the inverter is biased through the " high-gain" region and back again. 

The noise characteristics are most strongly affected by the body voltages on FETs 
Ml, M3, and M5. The diamonds on the left axis of Figure 10b 5 lib, 12b, and 13b indicate the 
"initial condition" body voltages used for FETs Ml, M3, M5, and M5, respectively. The actual 
(from dynamic simulation) and estimated maximum values (from the static noise analysis tool) for 
these body voltages are compared in Table VII. 



TABLE VII 


Transistor 


Actual Body 
voltage (state 5) 


Estimated max body 
voltage (state 5) 


M3 for case 2 


0.05 V • 


0.255 V 


M5 for case 3 


0:720 V . . s : '/ 


0.787 V 


M5 for case 4 


0.11 V 


0.255 V " 


Ml for case 1 ' .. 


0.'95 V 


0.97 V s 



In Case 1, transistor Ml 3 charges the internal node between transistors Ml and M2 
to V DD every precharge phase. This configuration would be common in bulk CMOS to mitigate 
charge-sharing noise but is potentially disastrous for SOI since it can allow the body of Ml to easily 
charge to a very high bias, nearly 1 V. When brought into state 150 by the noise condition, Ml 
draws bipolar leakage current. The high body voltage also reduces the threshold voltage, making 
the FET more sensitive to the coupling noise on the gate. 

Case 3 is similar, but the fact that net E is under steady switching means that both 
states 120 and 140 are accessible for FET M5. State 2 accessibility limits the maximum body 
voltage possible to only about 0.75 V. In this case, the body voltage is not high enough to result in 
significant parasitic bipolar leakage but the threshold voltage is still reduced enough to result in 
functional failure. 
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In Case 2, transistor Ml 4 discharges the internal node between transistors M3 and 
M4 every cycle. From the point of view of the modified accessibility analysis, this is equivalent to 
Case 4, in which both E and F are active. The conservative upper bound overestimates the body 
voltage. The main reason for this discrepancy is the neglect of state 160 which is exerting a 
5 downward pressure on the body voltage which is not (in the interests of conservatism) considered 
in the modified accessibility analysis. 

The foregoing merely illustrates the principles of the invention. Various 
modifications and alterations to the described embodiments will be apparent to those skilled in the 
art in view of the teachings herein. Gate leakage is anemerging influence on the body voltage and 
10 could easily be incorporated in to the estimation technique. It will thus be appreciated that those 
n skilled in the art will be able to devise numerous systems and methods which, although not explicitly 
l^j shown or described herein, embody the principles of the invention and are thus within the spirit and . t 
scope of the invention. ' : !. .•' "km.-. ■ 
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